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ABSTRACT 

Context. Observations of the CO J(l-O) 115 GHz and J(2-l) 230 GHz lines in comet C/1995 01 (Hale-Bopp) were 
performed with the IRAM Plateau de Bure interferometer on 11 March, 1997. The observations were conducted in both 
single-dish (ON-OFF) and interferometric modes with 0.13 km s~^ spectral resolution. Images of CO emission with 
1.7 to 3^^ angular resolution were obtained. 

Aims. The ON-OFF and interferometric spectra show a velocity shift with sinusoidal time variations related to the 
Hale-Bopp nucleus rotation of 11.35 h. The peak position of the CO images moves perpendicularly to the spin axis 
direction in the plane of the sky. This suggests the presence of a CO jet, which is active night and day at about the same 
extent, and is spiralling with nucleus rotation. The high quality of the data allows us to constrain the characteristics 
of this CO jet. 

Methods. We have developed a 3-D model to interpret the temporal evolution of CO spectra and maps. The CO coma 
is represented as the combination of an isotropic distribution and a spiralling gas jet, both of nucleus origin. 
Results. Spectra and visibilities (the direct output of interferometric data) analysis shows that the CO jet comprises ^ 
40% the total CO production and is located at a latitude ~ 20° North on the nucleus surface. Our inability to reproduce 
all observational characteristics shows that the real structure of the CO coma is more complex than assumed, especially 
in the first thousand kilometres from the nucleus. The presence of another moving CO structure, faint but compact 
and possibly created by an outburst, is identified. 

Key words. Comets: individual: C/1995 01 (Hale-Bopp) - Radio lines: solar system - Techniques: interferometric 



1. Introduction 

Millimetre spectroscopy provided many insights into the 
composition and physical properties of cometary atmo- 
spheres. Many cometary parent molecules issued from the 
nucleus were identified with this technique. The high spec- 
tral resolution capabilities, and the possibility to observe 
several rotational lines belonging to the same molecule, 
allowed us to retrieve unique information concerning the 
velocity and temperature of the expanding coma, and 
the anisotropy of gas production at the nucleus surface. 
Because of the low, diffraction-limited spatial resolution 
provided by radio dishes at millimetric wavelengths, at best 
typically 10^', studies concerning the spatial distribution of 
parent molecules in the coma were sparse. 

The exceptional brightness of comet C/1995 01 (Hale- 
Bopp) near its perihelion on 1 April, 1997, motivated a 
wealth of innovative cometary observations. Among them, 
interferometric imaging of rotational transitions of par- 
ent molecules was successfully attempted. The Berkeley- 
Illinois-Maryland Association (BIMA) array mapped comet 
Hale-Bopp in HON J(l-O) and OS J(2-l) with 9" angu- 



lar resolution. Spatial asyn imetries su ^ gestin^ the presence 
of gas jets were detected ([Veal et al.L l2000l : I Wright et al.L 
'19981 IWoodnev et~all . 12002). Constraints on the photodis- 
sociative scalelengths of HON and CS were obtained from 
the radial extent of their radio emissions (jSnvder et al.L 
.2001 ). Using the O wens Valley Radio Observatory (OVRO) 
millimetre array, I Blake et alJ (p^999) obtained maps of 
HON, DON, HNC and HDQ at 2-4'' spatial resolutions over 
2-3 h integration time. The presence of arc-like structures 
offset from the nucleus is reported for all species but HON, 
and interpreted in terms of jets of icy particles releasing 
unalterated gas contrasting with that outgassed from the 
nucleus. 



Interferometric observations of rotational lines in comet 
Hale-Bopp were also made with the Plateau de Bure in- 
terferometer (PdBI) of the Institut de Radio Astronomic 
Millimetrique (IRAM) at 1-3'' resolution. A short and 
preliminary account of these ob serv ations was g i ven i n 
IWink et all (|l999f ) JDesDoisl (|l999h and lHenrv et afl (|2002f ). 
Millimetre lines of CO, HON, CS, HNC, CH .sOH, H2S , 
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SO, H 2CO were mapped (|Wink et al.Lll999l : [BQissier et all 
Bockelee-Morvan, e-mail: l2007t ) . At the same time, continuum maps of the dust and 

nucle us thermal emissions were obtained ([AltenhoflP et al.L 
fT999h . The PdBI was also used in single-dish mode to 
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detect new cometary mole cules (iBockelee-Morvan et al.L 
120001 : ICrovisier et all l2004a^,1oV 

We present here observations of the CO J (1-0) (115 
GHz) and J(2-l) (230 GHz) lines performed on 11 March, 
1997 at the Plateau de Bure interferometer. Among the 
~20 molecules identified in comet Hale-Bopp, and more 
generally in cometary atmospheres, carbon monoxyde CO 
is of particular interest: 

1. This species is the main agent of distant cometary ac- 
tivity, as first evidenced in comet 29P/Schwassmann- 
Wac hmann 1 (Sen av fc: JewittI , 1 1994 ICrovisier et al., 
Il995t) . This was later confirmed in comet Hale-Bopp 
from its long-term monitoring, which showed the change 
from a CO-dominated to an H2 0-domina ted activity at 
heliocentric distances rh ~ 3-4 AU (Bive r et"aL| , |1997| , 
ll999aD . 

2. In comets within 3 AU from the Sun, CO is, most 
often, the second major gaseous component of the 
coma after water. CO production rates relative to wa- 
ter are highly variable from come t to comet, rang- 
ing; from less than 1 % to ^ 20 % ([Irvine et al.L l2000l : 
iBockelee-Morvan et al.L |2004 for a review) . That mea- 
sured in comet Hale-Bopp near perihelion is among 
the highest eve r observed in comets: ~ 20% (e.g., 
[Bockelee-Morvan et al.L |2000L [pTSanti et all . l200lf ). 

3. There is much debate on CO production mechanisms. 
Because CO has a low sublimation temperature, the nu- 
cleus surface is certainly deprived of CO ice. Therefore, 
CO should outgass at some depth inside the nu- 
cleus, possibly from pure CO ice sublimation and/or 
from amorphous water ice when crystallizing and re- 
leasing trapped molecules (e.g., lEnzian et al.L 1 19981 : 
ICapria et al., 2000, 2002), if indeed pre-cometary ices 
condensed i n amorphous fo rm, which is somewhat de- 
bated (Mousis et al.L l2000l ). Another mechanism pro- 
posed to explain the CO production of comet Hale- 
Bopp near perihelion is the release of CO trapped 
i n crystalline water ice during water ice sublimation 
([Capria et al.Ll2000l . [2QQ2[ ). Comparing the CO jets mor- 
phology, as the nucleus rotates, to those of less volatile 
species or dust might provide clues to the origin of CO. 

4. There are several observational evidences that a sig- 
nificant part of the CO observed in cometary atmo- 
spheres could be produced by a distributed source. 
From in situ measureme nts of the local CO de nsity 
in IP/Halley with Giotto. [Eberhardt et all ([1987D con- 
cluded that only about 1/3 of the CO originated from 
the nucleus. The spatial distribution of CO molecules 
deduced from i nfrared long-s li t obs e rvatio ns of comet 
Hale-Bopp led iDiSanti et all (Il999l . l200lD to suggest 
that one-half of the CO was released by a distributed 
source when comet Hale-Bopp was within 1.5 AU from 
the Sun. The spatial resolution of the CO maps ob- 
tained at PdBI, approximately 1000 km to 1700 km 
radius on the comet depending of the line observed, 
is below the estimated radial extension of the CO 
distributed source of ^ 5 x 10^ km ([PiSa nti et al., 
I2001L iBrooke et al.L 120031 ). Therefore, an important as- 
pect of the study of the CO PdBI interferometric 
data is that independent information concerning the 
existence of a CO distributed source in Hale-Bopp 
coma can be possibly obtained. The study of the ra- 
dial distribution of the CO molecules is presented in 



a separate paper ([Bockelee-Morvan fc BoissierL I2009D . 
The brightness distribution of both 115 GHz and 
230 GHz lines can be fully explained by pure 
nuclear CO production, provided that opacity ef- 
fects and temperature variations in the co m a are 
taken into account (iBipckelee-Morvan et al.L I2005L 
iBockelee-Morvan fc BoissierL l2009l ). 

The observations of the CO J(l-O) and J(2-l) lines, 
performed both in single-dish and interferometric modes, 
are presented in Sect. O They show evidence for a rotat- 
ing CO jet. A model simulating a CO spiralling jet when 
the nucleus is rotating is developed in Sect. [3l It allowed 
us to compute synthetic spectra, visibilities and interfero- 
metric maps as a function of time, to be compared to the 
observations. Observations are analysed in Sect HI and the 
model free parameters are retrieved. A discussion is given 
in Sect. O 

2. Observations 

2.1. Description 

Comet C/1995 01 (Hale-Bopp) was observed from 6 March 
to 22 March, 1997 with the Plateau de Bure interferometer 
of IRAM, located in the French Alps. The observations of 
the J(2-l) (230.538 GHz) and J(l-O) (115.271 GHz) CO 
rotational transitions were carried out on 11 March, from 
4h to 15h UT. On this day, comet Hale-Bopp was at the 
geocentric distance A = 1.368 AU and heliocentric distance 
Vh = 0.989 AU. The weather conditions were good to ex- 
cellent and the atmospheric seeing was ~0.4 for both 1.3 
and 3 mm receivers. 

The comet was tracked using orbital elements provided 
by D.K. Yeomans (JPL, solution 55). The ephemeris was 
computed by P. Rocher (IMCCE, Observatoire de Paris) 
with a program which takes into account planetary per- 
turbations. The first interferometric maps obtained on 9 
March (HCN J(l-O) line) revealed that both continuum 
and molecular peak intensities were offset by about 5- 
6'' North in declination (Dec) from the nucleus position 
provided by the ephemeris. Observations on March 11 were 
made with the ephemeris corrected by 6" North in Dec. 

The PdBI was used in the compact configuration CI 
(see Fig. [1]) with five 15-m antennas providing 10 baselines 
(the spacing between two antennas) ranging from ~20 to 
~150 m. In 1997, the PdBI comprised a fiexible spectral 
correlator made of six independent units, providing corre- 
lated spectra with 64 to 256 channels spaced by 0.039 MHz 
to 2.5 MHz. We used 256 channels of 78 kHz separation 
for the observations of the 230 GHz line, and 256 channels 
of 39 kHz separation for those of the 115 GHz line. The 
four other u nits were used for the continuum observations 
presented in lAltenhoff et al.l ([19991 ). The effective spectral 
resolution is a factor of 1.3 broader than the channel spac- 
ing, and corresponds to ~ 0.13 km s~^ for both lines. 

The observing cycle was : pointing, focusing, 4 min 
of cross-correlation on the calibrators (2200+420 BL Lac, 
MWC349 and 3C373), 2 min of autocorrelation, and 51 
one minute scans of cross-correlation on the comet inter- 
laced with scans on the phase calibrator (2200+420 BL 
Lac) observed every 20 min. The cycle was completed by 
another 2 min of autocorrelation on Hale-Bopp. For the au- 
tocorrelation observations (in this mode, the five antennas 
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Fig. 1. The Plateau de Bure Interferometer in CI configu- 
ration 
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Fig. 2. li-y-coverage at the PdBI on 11 March, 1997 with 
the CI configuration. For each arc of ehipse, the white cir- 
cle represents the ii'U-point at the beginning of the observa- 
tions. The time evolution of the li^'-points loci is counter- 
clockwise. 



behave as five independent single-dish telescopes), we used 
position-switching (ON-OFF) with 5' offset to cancel the 
sky background. The spectra of the five antennas were then 
co-added. Hereafter, these autocorrelation observations will 
be denominated as ON-OFF observations. 

The amplitude and phase calibrator was 2200+420. 
MWC349 was used to determine the flux density of 
2200+420. Bandpass calibration was made on 3C273. 
Because of less accuracy in phase calibration after 12.5 
h UT, only inter ferometric data acquired before 12.5 h 
UT were considered. Calibration was done with the IRAM 



CLIC software, and the data hence derived were stored in 
'^v-tables. Reduction and cleaning of the maps were per- 
formed with the MAPPING/GILDAS softwar^]. 

Concerning ON-OFF spectra, antenna temperatures 
were converted into main beam brightness temperatures 
TmB through TmB = ^A^eff/^eff with beam efficiencies B^f^ 
of 0.83 (115 GHz) and 0.58 (230 GHz), and forward efficien- 
cies Feff of 0.93 (115 GHz) and 0.89 (230 GHz). Flux density 
per beam {S in Jy) is then related to antenna temperature 
through S/Tl = {2knmB/X^)Feif/B,ff = IQMF.^/B.s, 
where ftmB is the main beam solid angle. 

For both ON-OFF and interferometric data, the uncer- 
tainties in ffux calibration are at most 10% and 15% for the 
J(l-O) and J(2-l) lines, respectively. The r.m.s in phase 
noise ranges from 10 to 27° at 230 GHz and from 4.6 to 
20° at 115 GHz, depending of the baseline. 

The cross-correlated spectra produce, for each spectral 
channel, interferometric maps with spatial resolutions cor- 
responding to the -i/i^coverage (Fig. [2j). When all cross- 
correlation data are considered, the full width at half max- 
imum (FWHM) of the synthesized beam is 2.00'' x 1.38'' 
with major axis at position angle pa = 99.56° at 230 GHz, 
and 3.58" x 2.57" with pa = 86.00° at 115 GHz. The 
FWHM of the primary beam of the antennas is 20.9" at 
230 GHz and 41.8" at 115 GHz. 

2.2. ON-OFF spectra 

ON-OFF spectra of the CO J(2-l) line are shown in 
Fig. [3^. The integration time is 2 minutes (on+off) for 
each spectrum. They show a feature moving from positive 
to negative velocities, and to positive velocities again, with 
respect to the nucleus velocity frame. In other words, a jet- 
like CO gas feature, which velocity vector with respect to 
Earth rotated during the course of the observations, is ob- 
served. This CO gas feature contributes up to 28% of the 
total line area. 

The synodic rotation period of comet Hale-Bopp in 
February-A pril 1997 has been dedu c ed from studies of the 
dust shells (ISarmecanic et al.L Il997l : iFarnham et al.L 1 19991 : 
lOrtiz fc Rodriguez) , 119991 ). Th e most accurate value P = 
11.31 ± 0.01 hr measured bv iFarnham et all (|l999f ) is in 
agreement with a slightl y larger sidereal ro t ation period 
P = 11.34 ± 0.02 hr (fLicandro et all . Il998l : iJorda et al.l . 
1 19991 ). Figure [4] plots the evolution of the line velocity 
shift (the spectrum first order momentum) with time. The 
points follow a sinusoidal curve which period corresponds 
to the comet's nucleus rotation (taken to be equal to 
11.35 hr throughout this paper). The sinusoid determined 
from a least-squares fit with fixed period P = 11.35 h 
has a mean level of vq = —0.05 ± 0.01 kms~^ and its 
amplitude is ^ = 0.29+0.03 kms"^ (Fig. Hj). The line 
area does not show significant variation with time, with 
a mean value of 4.22 + 0.03 Kkms~^ in main beam 
brightness temperature scale T^b (i-e., 82.8 Jy km s~^ in 
flux density scale). Fluctuations of ~10% at most are ob- 
served (with a standard deviation of 5%), which are not 
correlated with the velocity shift variations. The veloc- 
ity shift curve obtained for the J(l-O) line is much more 
noisy (error bars ~0.2 kms~^ on the individual spectra), 
but is similar to the J(2-l) velocity shift curve: a least- 
squares sinusoid fit with P fixed to 11.35 h leads to .4 = 

^ http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 3. CO J(2-l) spectra obtained in ON-OFF and interferometric modes as a function of time, given in UT hours 
on 11 March, 1997. The velocity scale is with respect to the comet rest velocity. Spectra have been shifted vertically 
according to observation time, but are scaled identically, a) ON-OFF spectra: the integration time is 2 min (on+off). 
b— d) Interferometric spectra (visibility amplitude as a function of spectral channel): averages of 25 to 27 scans of 1 min. 
Only data for the 3 shortest baselines are shown. 



0.39±0.16 kms-\ vo = -0.05±0.05 kms-^ (iHenrvl . 120031 ). 
Adding all spectra, the mean velocity shift of the J(l-O) 
line is —0.09 ± 0.05 kms~^, in agreement with that of the 
J(2-l) line (-0.083 ± 0.007 kms"^). The line area of the 
J(l-O) line is 0.552 ± 0.022 Kkms"^ in the T^b scale (i.e., 
10.8 Jy km s-^). 

From the spin axis orientation and the equatorial coor- 
dinates of the comet, it is possible to derive the angle 6^ 
(aspect angle) between the spin axis and the line of sight, 
and the North pole position angle pa^j , defined from North 
to East. The different spin orientations published in the 
literature are listed in Table [TJ 

Ad opti ng the spin orientation derived by iJorda et al.l 
(|l999l ) and I Schleicher et al.l (|2004 ), the comet spin axis was 
then only 20° far from the plane of the sky. In such a geo- 
metrical configuration, a polar CO jet would lead to an al- 
most constant velocity shift. A jet close to the equator can 
explain a velocity shift following a sinusoid centred around 
vq ^ km s~^. Both this sinusoidal curve and the constant 



CO line area show that the amount of CO gas released in 
this jet did not vary during nucleus rotation. Given the 
Sun direction (phase angle of 46°, pa = 160°), this near- 
equatorial CO jet was active night and day at about the 
same extent. 

From the line areas of the J(l-O) and J(2-l) ON-OFF 
profiles, we derive a CO production rate Qco = 2.1 xlO^^ 
s~^. Here, we have assumed a Haser parent molecule dis- 
tribution for CO, and run our excitation model (Sect. [3]) 
with a kinetic temperature T of 120 K which agrees with 
temperature determinations pertaining to the 10 000-20 000 
km (radius) coma region sampled by the primary beam of 
PdBI (|Biver et al.L Il999al : [DiSanti et alj, l200lf ). Using an 
extended production for CO consistent with the IR ob- 
servations does not significantly affect the inferred Qco- 
iDiSa nti et al.* ('2001') inferred a total CO production rate 
(nuclear+distributed) fully consistent with our value. In 
the following sections, we will assume T = 120 K and a 
total Qco of 2 xlO^o s-^ 
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Time on 11 March, 1997 [UT] 

Fig. 4. Time evolution of the velocity shift of CO J(2-l) 
ON-OFF spectra shown in Fig.[3K. The plotted curve is the 
least-squares sinusoid fit to the data. It has a fixed period 
of 11.35 h, an amplitude of 0.29±0.03kms~^, and a velocity 
centre vq = — 0.05±0.01 kms~^ (dotted line). 

Table 1. Spin axis orientation. Columns 1 and 2 are the 
equatorial coordinates found in the literature. Columns 3 
and 4 are the corresponding position angle (pa^) and aspect 
angle {0^) on 11 March, 1997. 
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2.3. Interferometric data 

We present in Fig.[5]the line integrated maps of CO J(2-l) 
and J (1-0). 8 hours of observations were averaged (from 
4.5 h UT to 12.5 h UT) and 25 velocity channels (12 on 
both sides the central channel corresponding to the nucleus 
velocity) were co-added. An asymmetrical shape, which 
is not aligned with the elliptical clean beam (the synthe- 
sized interferometer beam), is observed and related to the 
anisotropy of the gas emission. 

In the line integrated interferometric map of CO J(2- 
1) (Fig. [5]), the position of the peak brightness ( Cm) is 
at RA = 22h 30m 38.02s and Dec = 40° 46' 3.1'' (with an 
astrometric precision of 0.07") in apparent geocentric co- 
ordinates given for 7.00 h UT. The peak position of the 
CO J(l-O) brightness (RA = 22h 29m 38.46s and Dec = 
40° 41' 10.1" at 4.00 h UT) is consistent with that of J(2- 
1), taking into account the comet motion from 4 to 7 h 
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Fig. 5. CO J(l-O) (top) and J(2-l) (bottom) line inte- 
grated maps observed on 11 March, 1997 (all data). 25 spec- 
tral channels have been averaged. The synthesized beam is 
in the lower left. The dashed square on the J(l-O) map cor- 
responds to the size of the J(2-l) map. For the CO J(l-O) 
line, contours are 0.037 Jy/beam, and the r.m.s. is 0.018 
Jy/beam. For the J(2-l) line, contours are 0.186 Jy/beam 
and the r.m.s. is 0.066 Jy/beam. Iso-contours are successive 
multiples of 10% of the maximum intensity, at 10 to 100 % 
of the maximum intensity. The solar direction is indicated. 
The original images have been shifted so that the maximum 
brightness peaks at the centre of the maps. In units of line 
area, the intensities at maximum brightness are 0.93±0.04 
and 4.65±0.15 Jy km s"^ for the J(l-O) and J(2-l) lines, 
respectively. 




Fig. 6. CO maps as a function of spectral channel on 11 March, 1997 (all data averaged). RA and Dec positions are 
with respect to the mean photometric centre Cm determined from the whole data set. The velocity (with respect to the 
comet rest velocity) of the spectral channels is indicated in the top left corner of the maps. The synthesized beam is in 
the lower left, and is coloured according to the velocity value. J(l-O) line (top): contour interval is 0.053 Jy/beam and 
the r.m.s. is 0.046 Jy/beam. J(2-l) line (bottom): contour interval is 0.283 Jy/beam and the r.m.s. is 0.15 Jy/beam. 
Contours correspond to multiples of 10% the peak flux density measured on channels at Doppler velocities of —0.71 and 
-0.81 km s-^ 
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Fig. 7. Individual maps of CO J {2-1) for data subsets of 1 h. Iso-contours are successive multiples of 10% of the maximum 
intensity, at 10 to 100 % of the maximum intensity. For each map labelled (i) on the top left corners, a crosse locates the 
mean photometric centre Cm determined from the whole data set. The arrow represents the direction of the individual 
photometric centre d with respect to Cm- Ci is evaluated by fitting a 2-D Gaussian of adjustable width. The beam shape 
is shown in the bottom right corner. 
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Fig. 8. Time evolution of the photometric centres. 
Photometric centres with i refering to the maps shown 
in Fig. [71 are given with their errorbars. The ellipse drawn 
is a least-squares fit to the data. It corresponds to a paral- 
lel of latitude 10° North of a sphere, with spin axis aspect 
angle = 80° and position angle = 211°. 



UT (see iBoissier et al.L 120071 ). The peak position of the 
continuum emission at 230 GHz observed simultaneously 
also almost coincides {^-2" offset) with the C O J(2-l) peak 
(jAltenhoff et al.L Il999l : iBoissier et aD . I2007D . Using orbital 
elements based on optical astrometric positions from April 
1996 to August 2005 (JPL solution 220), the offset between 
the CO peak and the ephemeris is +2.9^' in Dec and +0.4'' 
in RA. So, positions of the CO and radio continuum bright- 
ness peaks differ by typically +3'' in Dec from optical as- 
trometric positions. A bright dusty jet was identified south- 
ward in th e optical images of comet Hale-Bo pp near perhe- 
lion (e .g., iJorda et^aH , Il999f ). As shown bv lBoissier et al.l 
(|2007f ). the optical astrometric positions were more af- 



fected by dusty jets than the radio positions. IBoissier et al.l 
(200^^ showed that the astrometric positions provided by 
the IRAM continuum radio maps provides an orbit which 
does not require the existence of non-gravitational forces 
acting on the Hale-Bopp nucleus, in contrast to those de- 
rived from only optical positions, thereby solving a con- 
tentious issue. In conclusion, there is no substantial offset 
between the nucleus position and the mean photometric 
centre of CO emission. 

The J(l-O) and J(2-l) spectral channel maps (Fig. [6]) 
were obtained with the same procedure. The peak bright- 
ness on the blue channels is stronger than that on the red 
ones. This indicates more emission toward the Earth and is 
related to the jet seen in ON-OFF spectra. Indeed, the in- 
terferometric observations covered only 2/3 of the nucleus 
rotation period, when the jet was, most of the time, fac- 
ing the Earth (Fig. [3|). The spectral maps show that the 
CO coma structure is complex. The interpretation of the 
brightness distribution on these maps is not straightfor- 
ward, since the signal is here averaged over the entire pe- 
riod of observation and the CO coma is rotating. The most 
central channels are sensitive to molecules expanding along 
directions close to the plane of the sky. They show coma 
structures towards North- West and South-East quadrants 
(roughly along a direction perpendicular to the projected 
rotation axis, see the 230 GHz maps in Fig.[6|). These struc- 
tures may trace the jet at the time it was near the plane of 
the sky. Channels at high negative velocities (-0.6 to -0.9 
km s~^) show a much brighter and strongly peaked inten- 
sity distribution because the jet is here facing the Earth. 

In order to investigate whether there is temporal evi- 
dence for the rotating jet in the images of the CO emission, 
we have combined the data into five parts of about 1 hour 
each. Resulting maps are presented in Fig. [71 Because of 
Earth rotation, the beam shape rotates with time from map 
to map and changes dimension (see next paragraph). This 
prevents a detailed study of the rotating jets directly from 
the maps and, as explained later, another approach will be 
used. However, an interesting feature is observed. From the 
observations averaged over the whole day, we have derived 
the mean photometric centre of CO emission. Cm- For each 
map z, we can also derive the photometric centre, Ci, and 
the vector J\ = CmCi, as shown in Fig. [71 The time evolu- 
tion of Ji is presented in Fig. [HI We observe that it moves 
counterclockwise (disregarding C5) along an ellipse which 
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long axis is perpendicular to the spin axis direction. Such a 
displacement is that expected in presence of a CO rotating 
jet. Provided Cm coincides with the nucleus position, Ji re- 
flects the jet direction on map i. For a spherical nucleus and 
constant jet activity, the Ci's locus should be then an ellipse 
which long axis position angle is perpendicular to the spin 
axis, the other characteristics of the ellipse (axis lengths, 
centroid) being related to the amount of CO gas inside the 
jet, as well as to its latitude on the nucleus surface. A least- 
squares fit of the photometric centres leads to an ellipse 
(Fig. [8]) that corresponds to a spin axis with position angle 
pa^ =211° and aspect angle 0^ =79°, in good agreement 
with most of the published values (Table [1]). The ellipse di- 
mensions and position inferred from the fit do not provide 
direct quantitative information on the jet relative strength 
and latitude because the nucleus position may be off by 
a fraction of an arcsec with respect to Cm- However, the 
significant displacement of the photometric centre during 
nucleus rotation excludes a high-latitude jet, in agreement 
with the conclusion obtained from the ON-OFF spectra. 
The small offset between Cm and the nucleus position (as 
determined from the peak of the continuum emission) is 
also consistent with a low-latitude jet. 

For a deeper study of the interferometric data, we have 
decided to work on complex visibilities in the i^t^plane. 
For readers not familiar with interferometry, let us explain 
briefly what this means and how maps are obtained. An 



interferometer measures the Fourier transform (FT) of the 
source brightness distribution on the sky. The complex visi- 
bilities V(ii, v) sample the FT at points (ix, v) in the Fourier 
plane, also called the i^-y-plane. These points are the pro- 
jections of the baselines onto the plane of the sky and de- 
fine the lii^coverage of the observations (Fig. [2]). As the 
Earth rotates, the locus of the points (u^v) produced by 
one baseline is an arc of ellipse. Therefore, the uv-ididms 
a = aA? + ^ changes with time (except if the source 
observed is circumpolar, because the locus is then a cir- 
cle). The longer we observe, the lo nger are these arcs, the 
larger is the uv-plane coverage (see [Thompson et al.L ll99lL 
chap. 4, §4.2). (In fact, the ^/i;- cover age produced by a cou- 
ple of antennas comprises two arcs of ellipse symmetrical 
with respect to the centre {u^v) = (0,0); this is because 
the source brightness distribution is a real function, so its 
FT verifies V(ii, v) = V(— li, —v).) In order to make a map, 
one has first to compute the inverse Fourier transform of 
the sampled signal. In a second time, this dirty map is de- 
convolved from the dirty beam, which is the FT of the uv- 
coverage. Because the ^/v-plane is not regularly covered, in- 
terpolations are made when performing the FT. In addition, 
when the i^-y-coverage is highly anisotropic, the dirty beam 
presents intense sidelobes, which might not be properly ac- 
counted for in the deconvolution step. This might result 
in the apparition of artefacts. The anisotropic -i/^^coverage 
also results in an elliptical clean beam. 
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Since the individual baselines have different directions 
and lengths, they probe different scales and regions of the 
coma. So, visibilities have to be studied for each baseline 
separately. We plot in Fig. [9] the time evolution of the vis- 
ibility amplitude V of the CO J(2-l) line with respect to 
the uv-T didius a. The visibilities have been integrated over 
velocity and have units of line area. Let us assume that the 
line is optically thin and its excitation does not vary within 
the field of view. For an isotropic coma described by a par- 
ent molecule distribution, the visibility curve would follow 
V(cr) oc cr~^, provided the photodissociation scale length is 
large compared to the f ield o f view , which is the case here 
(iBockelee-Morvan & B oissierl . I2QQ9[ ). We observe in Fig. [9] 
some modulations with respect to the mean evolution (in 
^-i.iSj ^YiQ,t are not due to noise. They cannot be due to 
variations of the activity of the comet since the area of the 
line, observed in ON-OFF mode, is roughly constant with 
time. Furthermore, modulations do not present the same 
behaviour from one baseline to another. 

We presented in Fig.[3)3-d spectra issued from the three 
shortest baselines of the interferometer. As on the ON-OFF 
spectra (Fig.[3K), we can see spectral features moving from 
red to blue velocities. Figure [TOl presents the time evolution 
of the interferometric velocity shifts. At least for the five 
shortest baselines, they can be fitted by sinusoids of period 
equal to the nucleus rotation period. We observe that these 
curves are not in phase. As the baselines do not have the 
same length, nor the same orientation, this phase difference 
may suggest that the CO jet is spiralling and detected at 
different times with the various baselines. A straight jet 
would produce velocity shift curves in phase. 

The next section presents a model with a rotating gas 
jet aimed at the interpretation of these data. 



3. Model 

The model is a static 3-dimensional model simulating a 
spiralling gas jet. It computes uv-tahles (i.e., visibilities) 
corresponding to the observational conditions of the data, 
that is with the same -i/t^ cover age. The time evolution of 
the coma is simulated by computing a serie of successive 
static i^i^tables. 

The CO coma is modelled as a combination of an 
isotropic outgassing and a gas jet defined by: 

— its half- width ^; 

— its latitude £ on the nucleus, assumed to be spherical; 

— the fraction fco of the CO released within the jet. 

The model parameters are summarized in Table [21 

The coordinate frame (Oxyz) used in the calculations 
has its origin at the comet nucleus with the z-axis along the 
line of sight opposite to the Earth, and the x and ?/-axes 
pointing East and North in the plane of the sky. tq is the 
nucleus radius. The jet direction at the nucleus surface Jo 
is {Oo^(j)o)t in the {Oxyz) frame at time t (usual spherical 
coordinates) and is moving with time due to the nucleus 
rotation. Let us define R^j the rotation matrix for a lapse 
At, so that we have (6>o,^o)t+A^ = Ru • (^o,0o)^. The jet 
direction at distance r = ro + Ve^pAt in the coma is Jr = 
R~^ ■ Jo, where 

'^exp is the gas expansion velocity. 
We assume a Haser-like parent molecule distribution for 
CO. Indeed, as commented in Sect. [U though infrared ob- 
servations suggest that part of the CO in Hale-Bopp coma is 
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originating from a distributed source (jPiSanti et al.L 120011 : 
iBrooke et al.L 120031 ) , t he present observations do not requir e 
CO to be extended ("Bockelee-Morvan & Boissier", I2009D . 
The local density at (r, 0, (j)) direction is then given by: 



ncoyr, 0, (j)) = -j^j^ exp 



{r - rp) 
^co 



where 



Q(r,^,(/)) =Qiso + 4^Qjeta(r,^, 



(1) 



(2) 



Qiso and Qjet are the total CO production rates due to the 
isotropic contribution and within the jet, respectively. For 
the sake of simplicity, the total CO production rate Qco = 
Qiso + Qjet is fixed, and taken to be equal to the 2 x 10^^ 
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s-^ (|Biver et alJ (|l999af ). see Sect.[22|). f^o = ^ is a free 
parameter. The function ^(r, 6, (j)) describes the jet pattern 
and is a normalized Gaussian (J Q(yt)dVt = 1) of half width 
^ centred on J^. Leo = Vexp/Pco is the photodissociative 
scalelength, where Pco is the CO photodissociation rate. 

The code uses = 47 (Oxy) grids, each of them 
sampling one channel of the spectrum centred on vi = 
{^^^ — i) Sv^ in the nucleus velocity frame, with Sv = 0.10 
kms-^ The (Oxy) grids are lOO^'xlOO'' long0. They are 
divided into 256 x 256 cells which dimensions are {6x^ Sy) 
{Sx = Sy = 0.39 

In the optically thin case, the brightness distribution in 
the plane of the sky [W m~^ sr~^], when selecting only 
molecules contributing to channel z, is given by: 



Table 2. Model parameters. 



Mi{x,y) 1 



(3) 



where v is the line frequency, and A^i is the Einstein co- 
efficient for spontaneous emission. 5x and 5y have units of 
radians. A is the geocentric distance. 

JVi{x^y) is the number of CO molecules in the upper 
state of the transition sampled by the cells with Doppler 
velocities contributing to channel i: 



nX+^ ny+^ n 10 Lco 

,y)= / / ncoPu^idxdydz, (4) 

A-4f Jy-^ ^-lOLco 



J\fi{x 



where Pu is the relative population of the upper level of 
the transition, which depends on the radial distance to nu- 
cleus, and nco is from Eq.[Tl H^, the function used to select 
velocities, is defined by: 



Ri{x,y,z) 



I 1 if^. 
I else^ 



{x,y,z) 
elsewhere 



Sv 1 
2 J 



(5) 



where Vz{x,y, z) is the gas velocity projected onto the line 
of sight. The gas velocity is radial, and its amplitude is a 
Gaussian centred on I'exp, which width is 2^hi(2)kT /m^ 
to account for thermal broadening, k is the Boltzman con- 
stant, T the kinetic temperature and m the CO molecular 
mass. 

Because the CO production rate is high in comet Hale- 
Bopp and a dense CO jet is present, optical depth effects 
need to be considered for the calculation of the CO J(2- 
1) brightness distribution (i.e., Fi{x^y)). They are not ex- 
pected to affect significantly the ON-OFF spectra, but 
could be significant for the inter ferometric signals. The re- 
sults presented in this paper were performed solving the full 
radiat ive transfer equation, as explained in iBoissier et al.l 
([2003), 

assuming the local velocity dispersion to be ther- 
mal. 

A synthetic 47-channels ON-OFF spectrum is obtained 
by the convolution of Fi with the antennas primary beam. 

For each channel i , the visibilities are defined by (see 
e.g.. [Thompson et al.L[l99ll . chap. 4, §4.1) : 



Vz(^) = ^ / A{s)F,{s) exp(-^o- • s)dQ, 
vdv J^^ c 



(6) 



this is much larger than the primary beam of the antennas, 
but this was necessary in order to have a good resolution in the 
Fourier plane. 



fixed parameters 



heliocentric distance 




0.989 AU 


geocentric distance 


A 


1.368 AU 


gas expansion velocity^^-* 




1.05 kms"^ 


total CO production rate*^^-* 


Qco 


2.0 X 10^° s-^ 


CO photodissociation rate*^^-* 


/3co 


7.50 X 10-^ 


nucleus rotation period 


P 


11.35 h 



free parameters 



fraction of CO in the coma outgassed in the jet fco 

jet latitude I 

jet half-width ^ 

rotation axis aspect angle Ouj 

rotation axis position angle pau 



iBiver et al.1 (|1999al) 
iHuebner et all ([199^ 1 



where cr is the baseline vector for two antennas, with coor- 
dinates {u^v) in the iiv-plane. s is a vector in the sky plane 
which coordinates are (x, y) in radian units. A is the power 
pattern of the antennas, and dVt is an element of solid an- 
gle on the sky. Vi(cr) is here in units of [W m~^ Hz~^] or 
janskys. Equation [6] can be approximated to : 



Vi(u,v) 



+CO z'+oo 



A{x,y)Fi{x,y) 



CO ^ — CO 



X exp( \ux + vy)) axay. 



(7) 



Visibilities are computed with a Fast Fourier Transform 
(FFT) algorithm. 

The population of the rotational levels {pu) is derived 
from an excitation model which takes into account colli- 
sions with H2O an d IR radiative pumping of the y (1-0) CO 
vibrational band ([Crovisier Le BourloU Il983l : IC rovisierL 
1987). This model provides populations Pu as a function 
of radial distance r, given a H2O density law with r. For 
simplicity, we assumed an isotropic H2O coma, and Pu only 
depending upon r. The collisional CO-H2O cross-sect ion is 
taken equal to dc = 2 x 10~^^ cm^ (jBiver et al.Lll999b[), and 
the H2O production rate is QU2O = 10^^ s~^ (Colo m et al.L 
1999). In the simulations presented in Sect.jH a kinetic tem- 
perature T = 120 K is used (see Sect. l2.2l for further discus- 
sion). The evolution of the population of the CO J = 2 and 
J = 1 levels is shown in Fig.[TTl The beam size of 20.9'' for 
CO J(2-l) corresponds to r ~ 10 000 km in the coma. Most 
CO molecules within the field of view are in local thermal 
equilibrium. 

Given the evolving coma, a full simulation of the ob- 
servations would require to compute for each one minute 
scan the visibilities corresponding to the current state of 
the coma and to their i^v-coverage. In order to limit the 
computer time, we modelled a whole nucleus revolution (P 
= 11.35 h) by 12 snapshots (see Fig. [12]). We tested the 
validity of the approach by verifing that calculations with 
an increased time sampling (namely 36 snapshots) provide 
similar results. Between 2 snapshots i and z + 1, the jet di- 
rection changed following (6>o,0o)i+i = Rui^o^Mi^ where 

is the rotation matrix for a span of P/12. Providing the 
initial jet longitude at time t corresponding to snapshot i=l 
is fixed, a composite i^t'-table can be computed with the uv- 
coverage of the observations. The numerical code computes 
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log distance from the nucleus [km] 

Fig. 13. Evolution of the amplitude of the velocity shift 
Fig. 11. Relative populations of the CO rotational levels curve A with O^j and ^, for fco = 66% and = 0°. The cross 
J = 1 and J = 2, as function of distance from the nucleus, indicates the value of ^ required to fit the observations {A 
A kinetic temperature of 120 K is assumed. = 0.29 ± 0.03 kms"^) for 0^ = 70-90° and fco = 66%. 




Fig. 12. Schematic view of comet Hale-Bopp nucleus on 11 
March, 1997, as seen from the Earth assuming pa^ =210° 
and Ouj =80° for the spin orientation. East is on the left. The 
latitudes are shown by steps of 10° and the arrow shows the 
rotation direction. The symbols represent the series of jet 
positions used to make a composite uv-tahle for an equa- 
torial jet. Plain (respectively dotted) symbols mean that 
the jet is on the visible (respectively hidden) side of the 
nucleus. 



twelve composite i^'U-tables, each of them corresponding to 
different initial jet longitudes spaced by 360° /12. The lon- 
gitude origin is chosen so that the sub-terrestrial point on 
the nucleus surface is at a longitude of 0°. For illustration. 
Fig. [12] shows the twelve jet positions for a jet at a longi- 
tude of 0° at the time of the first snapshot. The model also 
computes synthetic ON-OFF spectra for each snapshot. 



4. Jet morphology analysis 

In this Section, the ON-OFF and interferometric velocity 
shift curves observed for CO J(2-l) (Figs HI and [TQ|). and 
the time evolution of the visibilities (Fig. [9]) will be anal- 
ysed with the model presented in the previous Section to 
constrain its free parameters. Because of limited signal-to- 
noise ratio, it was not possible to make this analysis for 
J(l-O) observations. For the spin axis orientation defined 
by its aspect angle O^j and position angle pa^j , we restricted 
our study to the mean values found in the literature (see 
Table [J): 0^ = 60 to 90°, and pa^ = 200 to 230°. The jet 
width ^ was tested from 1° to 90°, the jet latitude £ from 
0° to 90° North, and f^o from 10% to 80%. 

4.1. ON-OFF velocity shift curve 

As discussed in Sect. 12. 2[ data shown in Fig.[4]are well fitted 
by a sinusoid with a period P = 11.35 h, an amplitude A = 
0.29 km s~^, and centred on = —0.05 km s~^. We can also 
define a phase to = 11. 75 ±0.12 UT, that corresponds to the 
time when the velocity shift is equal to vq on the increasing 
side of the curve. We used these three parameters {A^ vq^ 
to) as criterions for selecting the models that could explain 
the observations. Note that the position angle pa^ of the 
spin axis has no influence on the velocity shift. 

The phase to is only ruled by the initial longitude of the 
jet. From to derived from the observations, the best initial 
longitude is 300° at 3 h 47 UT. 

At first order, A is governed by fco and ^ : it increases 
when fco increases or ^ decreases. This behaviour is easily 
explained. If fco increases - all other parameters remaining 
unchanged - there is more signal coming from the jet falling 
into the same number of spectral channels. Then, the ve- 
locity shift increases, so does the amplitude of the velocity 
shift curve. In a similar way, when ^ decreases - keeping 
fco constant - an equal amount of signal coming from the 
jet is falling into a larger number of spectral channels. This 
results in reducing the velocity shift and also A. This is il- 
lustrated in Figs [13] and [TH which show the evolution of A 
with fco and ^. Moreover, A is not sensitive to the aspect 
angle except for small ^'s, as observed in Figs [13] and 
[T4l Furthermore, we note that the jet latitude i has little 
influence upon the amplitude, within the limits where we 
tested it {\£\ < 45°) (Fig. [T5|). To conclude, the observed 
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Fig. 14. Evolution of the amplitude of the velocity shift 
curve A with O^j and /co, for ^ = 20° and = 0°. The cross 
indicates the value of fco required to fit the observations 
{A = 0.29 ± 0.03 kms-^) for (9^; = 70-90° and ^ = 20°. 
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Fig. 15. Evolution of A (top) and vq (bottom) with ^ 
for several fco and £. Simulations have been done with 
0^ = 80° and pau; = 210°. Observed values are A = 0.29 
± 0.03 kms-^ and vq = -0.05 ± 0.01 kms"^ ) 



amplitude of 0.29 km s~^ can be fitted by many (/co,^/) 
combinations. 

The mean velocity vq of the simulated curves depends 
mainly on the jet latitude £. An equatorial jet always 
produces a curve centred on Okms~^, whatever Ouj. For 
0^ < 90°, a jet with a northern (resp. southern) latitude 
produces a curve centred on a negative (resp. positive) ve- 
locity. Again, this behaviour is understood. With O^j < 90°, 
the North pole is pointing towards the Earth (see Fig. [T2|). 
As a result, a northern jet is more often directed towards the 
Earth than a southern one. Opposite effects are obtained 
for 0^ > 90°, while 6^ = 90° (rotation axis in the plane of 



Table 3. Selected sets of parameters O^j, and fco re- 
producing the velocity shift curve of the ON-OFF observa- 
tions. 



set 




i 




fco 


(1) 


70° 


10° 


19.9° 


35.7% 


(2) 


80° 


10° 


11.9° 


33.6% 


(3) 


80° 


20° 


18.3° 


35.5% 


(4) 


80° 


30° 


33.0° 


47.8% 



the sky) always produces a curve centred on Okms~^, irre- 
spective of the sign of the latitude. The more the rotation 
axis is far from the plane of the sky, the more the veloc- 
ity shift curve is shifted. Furthermore, the velocity shift 
curve is all the more shifted as ^ is small and fco is large 
(Fig. [15]). This study leads us to the conclusion that, here 
again, many combinations (/co,^/) are able to reproduce 
the observed '^q. However, it shows that only a northern jet 
will be able to fit the observations. 

For a given iO^^/Jco) parameter set, it is possible to 
find the jet width ^ that is able to reproduce the ampli- 
tude A = 0.29 kms~^. Figure [16] shows in dotted curves the 
locus of the couples (/co,^) that reproduce the right ampli- 
tude A for several fixed (0^/) values. The same method is 
employed to determine the couples (/co,^) that reproduce 
the right velocity centre vq = —0.05 kms~^ (dashed curves 
in Fig. [16]). For each set (0^/)^ the intersection of the dot- 
ted and the dashed curves gives the only couple (/co,^) 
that reproduces A and vq. We made these computations 
for latitudes between 0° and 45°, and for 0^ = 70° and 
80°. Computations for O^j = 90° were useless because they 
provide vq = kms~^. 

The combinations {O^Jco^^/) selected by this study 
are summarized in Table [3] We note that the jet strength 
fco is typically between 35 and 50%. The jet is located on 
the northern hemisphere at a latitude between 0° and 45° 
(both excluded) for 0^ = 80°, and between 0° and 20° (both 
excluded) for 0^ = 70°. 

The velocity shift curve obtained for parameter set (3) 
is shown in Fig. [TO] together with the observed curve. We 
note the almost perfect match between model and observa- 
tions. The corresponding synthetic line profiles are shown 
in Fig. [TT] At any time, the jet contributes to both blue 
and red channels, with varying relative contributions, due 
to its spiral shape. 

4.2. Interferometric velocity shift curves 

In this Section, we study the velocity shift curves observed 
for the individual baselines. Figure [10] shows model re- 
sults with the parameter set (3) {6uj = 80°, pa^ = 210°, 
^ = 18.3°, £ = 20° and fco = 35.5%; Table[3]). Other param- 
eter sets of Table [3]give similar curves. Modelled curves are 
periodic functions, with a period equal to P. They mimic 
sinusoidal curves, though significant deviations from a si- 
nusoid are observed. This is because line shifts measured on 
V(ii, v) spectra are (i^, v) dependent: stronger jet contrast 
appears in specific regions due to spatial filtering. Then, 
due to the combination of Earth and jet rotation, regions 
with more or less jet contrast are sampled by the individual 
baselines. Simulations show that these curves evolve toward 
a true sinusoid when the jet width ^ is increasing {fco kept 
constant), due to smaller jet contrast. These curves change 
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when varying the jet parameters in the same way than does 
the ON-OFF curve. Changing the spin axis parameters by 
±10° does not affect much the curves. 

The modelled velocity shift curves for the different 
baselines are not in phase, with to (defining the phase, 
see Sect 14. ip increasing with decreasing baseline length 
(Fig. [18]). We expect a phase offset due to the spiral shape 
of the jet. Indeed, with respect to long baselines, short base- 
lines probe molecules in more distant regions of the spiral. 
Hence, they sample molecules released in average at ear- 
lier times. In addition, baselines of different length (even 
if they are parallel) record the maximum signal from the 
jet at different times due to the curvature of the jet. This 
can be understood from Fig. [191 which plots the amplitude 
of the visibility as a function of the orientation and length 
of the baselines for a simple geometry (rotation axis along 
the line of sight and equatorial jet) and at a given time. 
The combination of both effects introduces a phase offset 
in the velocity shift curves. The delay between two base- 
lines in the velocity shift curves represents the elapsed time 
between the jet detection by one baseline and its detection 
by the following one. Note that, given the large curvature 
of the spiral (molecules travel r ~ 10"^ km when the nucleus 
rotates by 90°), only its innermost part contributes signifi- 
cantly to the detected signal. The dashed curve in Fig. [T8l 
shows the evolution of to with 'U'^-radius, assuming that to 
varies linearly with a / Vexp • it follows approximately the to 
curve computed by the model (plain curve). 

The comparison between modelled and observed phases 
to and amplitudes A shows that there is relatively good 
agreement for some baselines and strong discrepancies for 
others (Fig. [T8|) . Good agreement for both to and A is ob- 
tained for baselines 3-4 and 2-4. Phase to is well repro- 
duced for baseline 1-3 (but not the amplitude). Strong dis- 
crepancy (by 3 h for to and a factor of almost 2 in A) is 
observed for baselines 1-4, 1-5 and possibly 4-5 (errorbars 
are large for this baseline). It is striking to note that good 
agreement is obtained for baselines with their field of view 
aligned along the spin vector, while discrepancies are ob- 
served for baselines with their field of view perpendicular 
to the spin vector. Clearly our model is too simple to repro- 
duce all observational characteristics. We did not find any 
simple explanation for these discrepancies. Velocity accel- 



eration inside the jet changes to evolution with baseline in 
the opposite way: it would result in a flatter increase with 
decreasing baseline length than obtained with a constant 
velocity. The presence of other CO evolving structures in 
Hale-Bopp coma is the likely explanation (Sect. [^3|) . 

4.3. Visibilities 

Figure [9] shows the time evolution of the visibilities V 
plotted as a function of the liv-radius a (as defined in 
Sect. 2.3, V refers to the amplitude of the visibilities in- 
tegrated over velocity). A least-squares fit to these data 
gives V(cr) (X cr-i i8±0 02^ compared to the varia- 

tion expected for a parent molecule distribution and an op- 
tically thin line (Bockelee-Morvan & Boissier, "2009). This 
trend can be explained by optical depth effects being more 
important for long baselines probing the inner coma. 

Modulations are observed around this fit: they trace 
variations of the brightness distribution sampled by the 
individual baselines as the baselines and jet are rotating. 
These modulations are characterized by their shape and 
their amplitude. The shape depends essentially upon the 
rotation axis position angle pa^j and the jet latitude i. For 
example, a high-latitude jet would result in strong modula- 
tions for baselines 3-4, 2-4 and 2-3, and no modulations for 
baselines 1-3 and 1-4 which scan regions along the equator 
(Fig. [2Q|) . The reverse is expected for a high-latitude jet 
and pa^ at 90° from the nominal pa^j of Hale-Bopp rota- 
tion axis. Qualitatively, the observed modulations exclude a 
high-latitude jet, as well as they exclude a rotation axis po- 
sition angle much different than the pa^j derived from visible 
observations. This confirms the conclusion obtained from 
the time evolution of the photometric centres (Sect. 12. 3p . 
sensitive to both the amplitude and phase of the visibilities. 

We compared the observed visibilities to those com- 
puted with the model with jet parameters (2), (3) and (4) 
given in Table [3] (those with O^j = 80° and pa^; = 210°). The 
lowest (reduced X%-n = 5.1 for A/" = 150 data points 
and n = 5 free parameters) was obtained for the set of pa- 
rameters (3) with ^ = 18.3°, i = 20° and f^o = 35.5%. 
Simulations with parameters (2) (respectively (4)) show 
larger (resp. lower) modulations than observed, and val- 
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Fig. 17. Synthetic CO 230 GHz ON-OFF line profile as a 
function of UT time on 11 March for parameter set (3) of 
Table [3] ((9^ = 80°, pa^ = 210°, ^ = 18.3°, i = 20° and 
fco = 35.5%). The total spectrum (jet+isotropic contribu- 
tion) is shown by thick lines. The grey spectrum shows the 
contribution of the CO jet. 



ues 60% (resp. 22%) larger than with parameters (3). Using 
jet parameters (3), we also made simulations with 0^ = 60 
and 70° and pa^j = 200, 220 and 230°. The was mini- 
mized for pa^j = 220° (Xat-u ~ 3.8), while has no sig- 
nificant influence on the visibilities. However, pa^j = 210° 
explains better the visibilities of the 3-4 baseline. Figure 
[2T1 shows the modelled visibilities with parameters (3) and 
pa^ = 210°. Looking to the shape of the modulations, there 
is an overall agreement between model and observations 
though, admittedly, the agreement is not perfect. A plot 
of observed versus modelled visibilities shows that the best 
agreement is for baselines 3-4, 1-5, 1-4, 1-3 and 4-5. The 
largest discrepancies are for baselines 2-4, 1-2, and 2-5. For 
most baselines, jet detection (traced by amplitude increase) 
occurs ~ 1 h earlier in the simulation than in the observa- 
tions. In contrast, in the velocity shift curves, the simulated 
jet is late with respect to the observed jet for most base- 
lines. This again shows that our model is too simple to 
explain satisfactorily the data. 

The visibilities obtained when optical depth effects are 
neglected show the same temporal behaviour (Fig. [22|) . 
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Fig. 18. Phase to (bottom) and amplitude A (top) of the 
velocity shift curves as a function of ^iv-radius. Dots with 
error bars are the data. The range of uv radius spanned 
by the baselines is shown with horizontal lines in the bot- 
tom figure. The plain curve shows model calculations with 
(9a; = 80°, pa^ = 210°, ^ = 18.3°, i = 20° and fco = 35.5%. 
The dotted curve in the bottom figure assumes that the 
elapsed time between jet detection by baseline at uv-ididms 
(Ji and jet detection by baseline at (72 (or in ON-OFF spec- 
trum at (72 = 15 m) is equal to ((72-(7i)/2'yexp, with Vexp = 
1.05 km s~^. The reference time is to measured on the ON- 
OFF velocity shift curve. The orientation of the baselines 
fringes during the course of the observations is shown at 
the bottom, together with the spin axis direction in the 
(RA, Dec) plane: e.g., the field of view of baseline 3-4 has 
its long dimension (reflecting the primary beam of the an- 
tenna) along the spin vector of the comet, and its short 
dimension (related to the baseline length) perpendicular to 
it. The Sun direction in the (RA, Dec) plane is also indi- 
cated. 

However, as expected, modulations are more contrasted in 
the optically thin case. 

The visibilities obtained with jet parameters (3) vary 
according to V{(j) cx cr"^*^^, which is consistent, in first 
order, with the observed variation (oc ^r"^-^^^^-^^). 

4.4. Maps 

Simulated maps as a function of time are compared to ob- 
served maps in Fig. [23l The shape of the observed CO 
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Fig. 19. Visibilities of the CO 230 GHz line (central chanel) 
as a function of baseline orientation in the uv plane and 
baseline length (a from 5 to 100 m by step of 5 m) at 
a given time. The rotation axis (P = 11.35 h) is along 
the line of sight {0^ = 0°). The jet is equatorial, has an 
aperture ^ = 30°, and its direction at the nucleus surface 
is at a position angle of 60° in the plane of the sky (i.e., 
Oq = 30°). Other parameters are Qco = 1 x 10^^ s~^, fco 
= 1, with assumed optically thin conditions. 
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Fig. 20. CO 230 GHz modelled visibilities for a high- 
latitude jet. 



coma is relatively well reproduced by the model. Some dif- 
ferences may be due to the presence of other CO coma 
features, as suggested by noticeable discrepancies at 8h50- 
9h50 UT. The time evolution of the photometric centre 
measured on the simulated data (jet parameters (3) with 
pa^ = 210°) is shown in Fig. [24l When the position of 
the photometric centre is refered to the mean photometric 
centre for the observing period, and is therefore directly 
comparable to the measurements, there is good agreement 
in the overall evolution. As could be expected, the relative 
(modelled— observed) positions generally differ, with dis- 
crepancies reaching -0.3'' for 4h30-5h30 and 7h20-8h20 
data (maps 1 and 3). However, the good overall agreement 
confirms a posteriori that the observed time evolution of 
the CO 230 GHz peak brightness position is related to the 
CO rotating coma. 
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Fig. 21. CO 230 GHz modelled visibilities for the jet and 

rotation axis parameters which reproduce at best the ob- 
servations (set (3) with pa^j = 210°). 
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Fig. 22. CO 230 GHz modelled visibilities for set (3) with 
pa^ = 210° assuming an optically thin line. 

Synthetic 230 GHz spectral maps directly comparable 
to the observations (Fig. [6]) were computed. Some observed 
basic features are well reproduced, as the very peaked and 
strong emission in the blue channels at high velocity and the 
elongated coma in East- West direction for the blue channels 
at velocities close to zero. However, the asymmetry of the 
spatial distribution observed in the red channels is not ob- 
tained in the model. As a matter of fact, the model, which 
parameters were constrained from the large scale ON-OFF 
observations, predicts that the inner parts of the jet probed 
during the course of the interferometric observations were 
most of the time projecting Earthward (i.e., with negative 
Doppler velocities). The discrepancies between modelled 
and observed spectral visibilities discussed previously are 
observed directly on the spectral maps. 

4.5. Evidence for a second moving structure 

The inability of our one-jet model to reproduce satisfac- 
torily the interferometric data is due to the presence of a 
second moving CO structure, possibly produced during an 
outburst. 

This moving structure is that seen North- West from the 
nucleus at positive velocities in the time averaged channel 
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Fig. 25. Time evolution of the CO 230 GHz photometric centre for individual velocity channels at v = -1.0, -0.5, 0, 
+0.5, +1.0 km s~^. The spectral data were smoothed to a spectral resolution of 0.5 km s~^. The time intervals labelled 1 
to 6 are given in Sect j4.5l . Note that they do not correspond to those used for Figs. [231 and [Ml The top figure shows the 
measurements, while the model (set of parameters (3) in Table [3j) is shown in the bottom figure. In order to illustrate 
how the photometric centres are affected by the main jet, we have connected points (coloured in blue) indicative of the 
motion of the CO main jet. The so-called "red jet" is shown by red symbols. 



maps (Fig. [6]). This structure, which is moving away from 
the observer, is oriented along the fringes of baselines 1- 
4, 1-5, 3-5 and 4-5, i.e., along those for which a strong 
discrepancy in the time variation of the velocity shift is ob- 
served. This structure is essentially detected by these base- 
lines : hence, the velocity shift measured for these baselines 
is smaller than expected (Fig. [TQ|). 

In order to study the time evolution of this stucture, 
the spectral data were smoothed to a spectral resolution of 
0.5 km s~^, with velocity channels centred at -1.0, -0.5, 0, 
+0.5, +1.0 km s~^ with respect to the comet rest velocity. 
Five consecutive time intervals (#1, #2, #4, #5, #6) of 
~ 50 min to lh30 long, were considered. The time inter- 
vals correspond to characteristic line shapes for the spec- 
trum recorded by baseline 3-4 (Fig. [3]): #1 (4h30-5h48 UT) 
symmetric shape ; #2 ( 5h48-7h07 UT) and #4 (7h07-9h08 
UT) asymmetric shape with strong blueshift ; #5 (9h08- 
lOhOO UT) symmetric shape; #6 (lOhOO UT-llh28) asym- 
metric shape with high redshift. A sixth period, labelled 
#3, merging intervals #2 and #4, was also considered. 

The 5x6 extracted data sets were analysed in the 
Fourier plane assuming that the brightness distribution is 
a 2-D Gaussian, with its intensity, width and offset with 
respect to the line-integrated, time- aver aged centre Cm as 
free parameters. Figure [25l shows the observed time evo- 
lution of the offset (RA(t), Dec(t)) for each channel. The 
same plot is shown using model (3) as input. 

For a non-rotating nucleus, the motion of a cloud of gas 
released during a short time interval from a nuclear source 



would be seen by a straight line in the (RA, Dec) quadrant 
and the velocity channel where its velocity vector projects. 
For a rotating nucleus, a permanent nuclear source near the 
equator will produce an expanding jet that will appear suc- 
cessively in the different channels as the nucleus rotates. In 
the case of comet Hale-Bopp, if the equatorial jet projects 
North- West in (RA, Dec) at positive velocities channels 
(V-^) at a given time, then it will successively project North- 
West in negative velocity channels (V-)^ South-East/VL, 
South-East/y^. As shown in Fig.[25l there is a satisfactory 
agreement between the model and the observations for the 
channels centred on -1.0, -0.5 and km/s. We have tenta- 
tively connected the points corresponding to the motion of 
the CO main jet. At time #1, the jet is Westward in the 
V = channel, at time #4 it is moving Eastward in the v 
= -0.5 km s~^ channel, at time #5 Eastward in the v = 
channel, and at time #6 Eastward in the v = +0.5 km s~^ 
channel. 

The motion of the red structure is clearly apparent on 
the channels at ^' = +0.5 and +1 km s~^. The direction of 
the motion suggests that it originates from a low latitude 
region at the nucleus surface. The longitude of this source 
is estimated to be within 90-150° westward from the source 
of the main jet. The weak contribution of the red structure 
to some baselines can be explained by spatial filtering, im- 
plying a compact structure. The CO 230 GHz fiux density 
in channels v = +0.5 and +1 km s~^ is four times smaller 
than in channel v = -0.5 km s~^ at the time the CO main 
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Fig. 23. CO 230 GHz Hale-Bopp maps as a function of 
time. Left column: observations; central column: simula- 
tions with the jet and rotation axis parameters which repro- 
duce at best the observations (set (3): d^j = 80°, pa^; = 210°, 
^ = 18.3°, i = 20° and f^o = 35.5%); right column: simu- 
lations for an isotropic coma. 

jet is contributing. This may explain why the red source 
does not contribute much to the ON-OFF spectra. 

5. Discussion 

5.1. Model assumptions 

The observations presented in this paper were interpreted 
with a simple geometric model of a CO rotating coma. 
Besides the assumed conical shape of the jet at ejection 
(see the discussion below), several other simplifying as- 
sumptions were made to limit the number of free model 
parameters. The outflow velocity was fixed and taken to be 
constant throughout the coma, where as some acceleratio n 
is expected by gas dynamics models (jCombi et al.l . [1999). 
Day-to-night asymmetries in velocity for the CO back- 
ground gas were not taken into account either, though they 
are clearly present. Indeed, the velocity cutoff and width of 
the blue wing of the CO velocity profiles are 10% higher 
than the corresponding values for the red wings (Fig. [3|), 
which suggests higher velocities towards the sunlit side of 
the nucleus given the Earth-comet- Sun geometry. In addi- 
tion, we assumed that CO molecules inside the jet expand 
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Fig. 24. Time evolution of the photometric centre, as deter- 
mined from fitting in the uv-plane. Measurements on Hale- 
Bopp CO 230 GHz observations are shown by black dots 
with error bars. Those from simulations with jet model (3) 
{pciuj = 210°) are plotted with red squares: empty squares 
show absolute positions; filled squares show photometric 
centres with respect to the mean photometric center over 
the period, which can be directly compared to the real data. 
The accuracy of the photometric centre position measured 
on the simulated data is also shown. 



at the same velocity than those in the background. We be- 
lieve that including slight variations in the flow velocity 
field would not change the main conclusions of this paper, 
though the characteristics of the jet could somewhat differ. 
For example, including this day-to-night asymmetry in ve- 
locity would shift the jet towards lower latitudes. Higher ve- 
locities inside the jet would require a proportionally higher 
jet contribution to the total CO production fco to fit the 
ON-OFF data. We also did not consider possible tempera- 
ture variations in the coma which have direct effects on CO 
excitation. A kinetic temperature higher in the jet than 
in the background gas would have resulted in a higher in- 
ferred fco' However, a 10 K difference in temperature would 
change fco by less than 10%. 

5.2. CO jet 

Our observations, interpreted with the help of a simple 
model, suggest the presence of a spiralling "jet" of CO, 
originating from the nucleus. The rotation of this jet is 
consistent with the rotation period and axis direction de- 
rived from most optical studies. The jet originates from 
a low-latitude (~ 20°) region of the nucleus, and com- 
prises a significant amount of the total CO production (~ 
40%) in a ~20-30° wide aperture. This is the first evi- 
dence of a CO spiralling coma around a comet nucleus. 
CO ro-vibrational line emissions observed in comet Hale- 
Bopp from January to May 1997 show strong East/West 
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asymmetries, with the position of the maximum brightness 
moving with time from up to ^2'' West to ^2" E ast of the 
continuum peak (DiSa nti et all [2QQll : iBrooke et aL. 2003). 
This is likely caused by the rotating CO structure observed 
in the radio ([Brooke et all I2QQ3I ) and shows that this CO 
structure was not a transient phenomenon related to some 
outburst. 

Our model, which includes a single jet, appears to be 
too simple to interpret all observation characteristics. Other 
"jets" are present, in particular another moving structure 
has been identified at positive velocities in the interferomet- 
ric data. Its full characterization is difficult from our data 
and beyond the scope of this paper. The CO main jet is cer- 
tainly more complex than a simple conical spiralling struc- 
ture. Gas structures created by inhomogeneities in gas pro- 
duction at the nucleus surface will unescapably be modified 
during their outfiow by the surrounding collisional environ- 
ment (the collisional sp here in Hale-bopp coma exceeded 
10^ km near perihelion. ICombi et allfTQQOl ). 

Gas dynamics models investigating the behavior of ex- 
panding dusty gas j ets have been developed (see the review 
of lCrifo et allf2005h . One or several active sources releas- 
ing water were considered, the rest of the nucleus being 
assumed to be weakly active. They show that, due to lat- 
eral expansion driven by pressure gradients, jets broaden 
with increasing distance to nucleus and interact with each 
other, to become eventually almost indiscernible in the 
outer coma. In addition, if there is solar modulation of the 
surface gas fiux, the jet morphology will change its appear- 
ance as the nucleus rotates. So far, these calculations may 
not be relevant to the CO jet seen in comet Hale-Bopp. 
First of all, it is very likely that CO diffuses from subsurface 
layers. Our observations suggest that, in ffrst order, there 
is no significant day-to-night asymmetry in the amount of 
CO gas released by the active source. Second, if the active 
source responsible for the CO jet does not produce similar 
strong enhancement in water production, pressure smooth- 
ing will not be acting much and jet broadening may rather 
result from molecular diffusion. In this case, it is possible 
that a nearly invariant CO density pattern co-rotates with 
the nucleus and remains relatively well preserved at large 
scales in the coma. However, it can be anticipated that the 
jet development will be still inffuenced by the surrounding 
environment. 

A good illustration is the following. iRodionov fc Crifol 
(|2QQ5f ) performed time-dependent 3-D gas dynamics cal- 
culations simulating rotationally-induced gas coma struc- 
tures. They used an arbitrarily Halley-like aspherical nu- 
cleus, scaled to Hale-Bopp nucleus size and homogeneous 
in composition. CO was assumed to diffuse uniformly from 
below the surface, while H2O sublimates according to so- 
lar illumination. The computed CO coma is asymmetric. 
The model produces faint H2O and CO spiralling struc- 
tures resulting from weak shocks induced by the surface 
topography. In other words, the CO outffow is largely in- 
ffuenced by the general and detailed properties of the ffow. 
The analysis of the CO Plateau de Bure observations using 
time-dependen t (multi-fluid) gas dynamic s calculations was 
performed bv .Boissier et al.! (|2QQ5l . I2QQ9D . They show that 
the observed time variations cannot be explained by the 
above mentioned shock structures and is due to a strong 
inhomogeneity in CO production from the nucleus surface 
or sub-surface. It is possible to hypothesize lengthly about 
the meaning of this CO overproduction. It can traduce in- 



homogeneities in CO content inside the nucleus but also lo- 
cal variations in mantle thickness, or in some properties of 
the nucleus material (e.g., dust/ice matrix structure, ther- 
mal conductivity). We leave this discussion to experts and 
encourag e them to perform nu merical simulations (see the 
review of lPrialnik et alll2QQ5[ ). 

5.3. Comparison with other studies 

Two questions arise: 

1) Is the CO jet we observed related (correlated) to jets 
observed for other molecules? 

2) Is the CO jet we observed related (correlated) to 
observed dust jets? 

In their interferometric observations of several molecu- 
lar species (HCN, H NC, DCN, HDO) with the OVRO array, 
iBlake et al.l (|1999[ ) observed that the molecular emissions 
peak at positions offset by a few arcsec from the contin- 
uum emission of the comet, and attributed these offsets 
to molecular jets. We do not observe such large offsets for 
CO when averaging the data over 70% of the rotation pe- 
riod. However, offsets reaching 1-1.5'' were observed for 
other molecules observed at the PdBI, interpreted to be the 
gaseous signa tures of the high la titude dust jet observed in 
the visible (Boissier et aLVSOOT*). The strong CO jet has no 
strong H2S, CS, SO and HCN counterparts though weak ro- 
tational modulations in the line sha pes of their radio em is- 
sio n are observed for so me of them ("Boiss ier et al.L [2QQ7I). 

iLederer fc Campinsl (2002) and Lede rer et all (|2009[ ) 
observed spiralling jets of OH, CN and C2. Several (up to 
five) active areas were necessary to reproduce their obser- 
vations from their Monte Carlo simulation. The strongest 
jets are coming from low-latitude regions (-22° to +20°). 
One of them could thus be associated with our CO jet. One 
low-latitude southern area is responsible for one-half of the 
OH coming from jet-source, but has a large opening angle 
(120°). So it is likely not related to the CO jet. Spatial 
profiles of H2O obtained from long-slit observations of in- 
fra red water lines only sh ow small East- West asymmetries 
(|Dello Russo et all l2000f ) . In summary, there is no strong 
evi dence for a stron g H2O jet associated to the CO jet. 

IWoodnev et al.l p002) mapped HCN in comet Hale- 
Bopp with the BIMA array. Their observations, which 
spanned over several days, were binned and averaged ac- 
cording to the phase of the comet rotation, in order to avoid 
smearing due to rotation. They found a jet morphology, 
possibly correlated with CN imaged in the visible. However, 
there is apparently no correlation of HCN with dust jets. 
The presence of a high-lati tude HCN jet is how ever sug- 
gested from the PdBI data (jBoissier et alll2007[ ). 

Dust jets also appear to be uncorrelated with our 
CO jet. Indeed, from vi sible imag ing at the time 
of the PdBI observations, iJorda etal.1 (|1999l ) observed 
a high latitude (64°) du s t jet . In the analysis of 
IVasundhara fc Chakrabortvl (| 19991 ). strong dust jets are 
present from latitudes +65°and -65°, in a ddition to weake r 
jets from low latitudes (see also Schleic her et all 12004 ). 
In February-May 1997, the high-latitude northern jet pro- 
duced repetitive sunward shells instead of a full spiral 
(Fig. [26|) . indicating that the source of this dust jet shut 
off during night time in contrast to the CO jet. 

In addition to the CO main jet, we have identified a 
second moving structure in red channels that we believe 
originates from a low latitude source. Interestingly, on 12.14 
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Fig. 26. Drawings of the inner coma of comet Hale-Bopp 
on 11 and 12 March, 1997 (N. Biver, 25.6 cm Newtonian 
telescope). 



March 1997 UT the dust coma presented a well defined shell 
at ~ 4C in a direction opposite to the repetitive sunward 
shells produced by the high-latitude sources (Fig. [26]). This 
shell was not present on 11.16 March, nor on 14.80 March, 
and can be attributed to an outburst initiated on 11.20 
March UT, i.e., at the time when the red structure began 
its expansion. Possibly, the red structure is related to this 
outburst. 

5.4. Plans for further studies 

ON-OFF and interferometric data were obtained at the 
Plateau de Bure interfer ometer for HCN, B[NC, CS, H2S , 
SO H 2CO and CH3OH (|Wink et all . Il999l : iBoissier et all . 
I2OO7D . Some lines (e.g., lines of HCN and CS) do show 
rotation-induced variations in their velocity shifts. The 
analysis of the Plateau de Bure observations is continuing 
and will be presented in forthcoming papers. 

The present study showed that radio observations can 
provide valuable information on the distribution of parent 
molecules in inner cometary atmospheres and its temporal 
evolution. In contrast with standard imaging techniques, 
radio observations are sensitive to radial velocities, i.e. they 



are sensitive to the gas distribution along the line of sight, 
whereas the other techniques are rather sensitive to the dis- 
tribution on the plane of the sky. They also probe different 
gas species. Radio observations and other techniques are 
therefore complementary. 

Radio interferometric imaging is a powerful tool for as- 
trometry. Our observations show that, in addition, rotat- 
ing comas can be detected from the motion of the centroid 
of molecular maps. This opens new perspectives because 
possibly useful constraints on the rotation properties of 
cometary nuclei will be obtained from such measurements. 

The analysis of our interferometric data was hampered 
by to the limited instantaneous li'U-coverage of the Plateau 
de Bure interferometer. The Atacama Large Millimeter and 
submillimeter Array (ALMA), with its 50 antennas, will be 
able to obtain images of molecular and continuum emissions 
with a short sampling time, high sensitivity and high an- 
gular resolution. It will provide a 3-D dynamical picture of 
inner cometary gaseous atmospheres, simultaneous images 
of the du st conia, and spatial informatio n on the gas tem- 
perature (iBiverl . [2QQ5I : iBockelee-Morvani I2OO8). Important 
breakthroughs concerning nucleus and coma processes can 
be expected. 
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